Nucleon dynamics in effective field theory of nuclear forces 
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Low energy nucleon dynamics in the effective field theory (EFT) of nuclear forces is investigated 
by using the formalism of the generalized quantum dynamics (GQD). This formalism is based on a 
generalized dynamical equation derived as the most general equation of motion consistent with the 
current concepts of quantum theory and allows one to extend quantum dynamics to the case where 
the fundamental interaction in a quantum system is nonlocal in time. The Schrodinger equation 
follows from this equation in the case where the interaction in a system is instantaneous. We 
show that the effective theory of nuclear forces gives rise to low energy nucleon dynamics which is 
governed by the generalized dynamical equation with a nonlocal-in-time interaction operator. This 
operator is shown to provide a natural parametrization of the interaction of nucleons and can be 
derived order by order by using methods of subtractive EFT's. We show that the use of the GQD 
for describing nucleon dynamics allows one to formulate the EFT of nuclear forces as a completely 
consistent theory with a well-defined equation of motion which does not require renormalization. 



I. INTRODUCTION 

Effective field theories (EFT's) are an important tool 
for computing physical quantities in theories with dis- 
parate energy scales [1-2]. Following the early work of 
Weinberg and others [3-7], the EFT approach has be- 
come very popular in nuclear physics [8,9]. To describe 
low energy processes involving nucleons and pions, all 
operators consistent with the symmetries of QCD are in- 
cluded in an effective Lagrangian. A fundamental diffi- 
culty is that such a Lagrangian yields graphs which are 
divergent, and gives rise to singular quantum mechanical 
potentials. To resolve this problem one has to use some 
renormalization procedure which regulates the integrals 
and subtracts the infinities in a systematic way. How- 
ever, as is well known, there are not any equations for 
renormalized amplitudes in subtractive EFT's. For this 
reason in some EFT's (see, for example, Refs.[6,10,ll]) 
finite cut-off regularization is used. Nevertheless, even 
when it is assumed that the cut-off has a physical mean- 
ing renormalization is required to render such theories 
consistent, and certain cut-off-dependent terms have to 
be absorbed into theories constants before determining 
them from empirical data. 

Another problem of the EFT approach is that one 
cannot parametrize the interaction of nucleons, by us- 
ing some Lagrangian or Hamiltonian. In fact, knowing 
the effective Lagrangian is not sufficient to compute re- 
sults for physical quantities. In addition, one needs to 
specify a way to make the physical predictions finite. 
The counterterms of renormalization involve additional 
unknown parameters. A renormalization scheme allow 
one to use experimentally determined parameters instead 
of the above unknown ones. Thus for computing phys- 
ical quantities one needs not only knowing the effective 
Lagrangian of the theory but also knowing some exper- 
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imentally determined parameters. The nucleon- nucleon 
(NN) interaction cannot also be parametrized by using 
the singular potential produced by the above effective 
Lagrangian: These potentials do not make sense with- 
out regularization and renormalization. On the other 
hand, if the EFT approach, as it is widely believed, is 
able to provide a fundamental description of the inter- 
actions of nucleons at low energies, one can hope that 
it will give rise to parametrization of these interactions 
by an interaction operators being so fundamental as the 
Coulomb potential that parametrizes the interaction of 
charged particles in low energy QED. In the quantum 
mechanics of particles interacting via the Coulomb po- 
tential, which is an example of an effective theory, one 
deals with well-defined interaction Hamiltonian and the 
Schrodinger equation governing the dynamics of the the- 
ory. This theory is internally consistent and provides an 
excellent description of atomic phenomena at low ener- 
gies. It is natural to expect that the EFT of nuclear 
forces is also able to provide such a fundamental descrip- 
tion of nuclear phenomena at low energies. In particular, 
one can expect that this theory will allow one not only to 
calculate scattering amplitudes but also to construct the 
evolution operator describing the dynamics of nucleon 
systems. In the present paper we show that recent devel- 
opments in quantum theory permit such a formulation of 
the effective theory of nuclear forces. 

The above problem of subtractive EFT's is the same 
that arises in any quantum field theory with UV diver- 
gences: Regularization and renormalization allow one to 
render the physical predictions finite, however, it is im- 
possible to construct a renormalized Hamiltonian acting 
on the Fock space, i.e. after renormalization the dynam- 
ics of the theory is not governed by the Schrodinger equa- 
tion. At the same time, in Ref.[12] it has been shown 
that the Schrodinger equation is not the most general 
dynamical equation consistent with the current concepts 
of quantum physics, and a more general equation of mo- 
tion has been derived as a consequence of the Feynman 
[13] and canonical approaches to quantum theory. Be- 
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ing equivalent to the Schrodinger equation in the case 
of instantaneous interactions , this generahzed dynami- 
cal equation permits the generalization to the case where 
the dynamics of a quantum system is generated by a 
nonlocal-in-time interaction. It has been shown [12] that 
a generalized quantum dynamics (GQD) developed in 
this way provides a new insight into the problem of UV 
divergences. 

It should be noted that in the Hamiltonian formalism 
the interaction generating the dynamics of an isolated 
system cannot be nonlocal in time in principle. In fact, 
such a nonlocality of the interaction should be associated 
with an energy dependence of the interaction Hamilto- 
nian, despite Hamiltonian itself is an operator represent- 
ing the total energy of the system. This conceptual prob- 
lem manifests itself in the fact that energy-dependent 
Hamiltonians are not Hermitean, and, as a result, the 
evolution of the system is not unitary. Such peculiari- 
ties are typical in definitions of renormalized interaction 
Hamiltonians (see, for example, Refs. [14,15]). The rea- 
son for this is quite obvious. The Schrodinger equation is 
local in time, and the interaction Hamiltonian describes 
an instantaneous interaction. Hence the introduction of 
a nonlocal-in-timc (energy dependent) interaction Hamil- 
tonian into the Schrodinger equation leads to an internal 
inconsistency of the theory. A remarkable feature of the 
generalized dynamical equation is that it provides the 
extension of the theory to the case of such interactions. 

Quantum mechanics is one of the basic ingredients of 
quantum field theory (EFT). For this reason one may 
expect that in nonrelativistic limit the Schrodinger equa- 
tion can be used as an equation of motion in the EFT of 
nuclear forces. However, this is not true, since within the 
framework of this theory the Schrodinger equation does 
not make a since without regularization and renormal- 
ization. Meanwhile, as has been shown in Ref.[12], only 
the generalized dynamical equation must be satisfied in 
general. This gives us the hope that this equation can 
provide a perfectly satisfactory description of low energy 
nucleon dynamics in the effective theory of nuclear forces. 
In the present paper we show that in leading order of the 
EFT approach low energy nucleon dynamics is governed 
by the generalized dynamical equation with a nonlocal- 
in-time interaction operator. Moreover, this dynamics is 
just the same as in the case of our model with a separa- 
ble interaction operator [12,16] that was used as a test 
model illustrating the possibility of going beyond Hamil- 
tonian dynamics provided by the GQD. As we show in 
the present paper, in the case of the quantum mechanics 
of nucleons at low energies, one has to deal with such 
a non-Hamiltonian dynamics. In Sec. Ill this is proved 
precisely in leading order of the EFT approach. We will 
show that the leading order contact component of the 
two-nucleon T-matrix obtained by summing bubble dia- 
grams, coincides with the T-matrix of the above model 
in the particular case where the form factor in the sepa- 
rable interaction operator is of the form (p{p) = 1. This 
T-matrix has the properties that are at variance with 



the ordinary requirements of quantum mechanics, and 
does not satisfy the Lippmann-Schwinger (LS) equation. 
At the same time, this T-matrix is the solution of the 
generalized interaction operator with a nonlocal-in-time 
interaction operator, and this operator parametrizes the 
contact component of the NN interaction in leading or- 
der. 

In Sec. IV we will investigate the features of low en- 
ergy nucleon dynamics in leading order of the EFT ap- 
proach. The possibility to extend our results to higher 
order will be investigated in Sec.V by using the example 
of the EFT of short-range forces developed by van Kolck 
[Q. We will show that in next-to- leading order, the dy- 
namics of this theory is also governed by the generalized 
dynamical equation with a nonlocal-in-time interaction 
operator. By using this operator as an example, we will 
demonstrate how one can construct order by order the 
interaction operator parametrizing the NN interaction in 
the effective field theory of nuclear forces. The new pos- 
sibilities that the GQD opens for practical calculations 
in the EFT approach will be discussed in Sec. VI. The 
advantages of this formalism become apparent in the case 
where numerical calculations are needed, i.e. in the case 
where the long-range components of the NN interactions 
are taken into account. It will be shown that for solving 
the evolution problem in this case the generalized dy- 
namical equation can be reduced to an integral equation 
that does not require regularization and rcnormalization, 
and is so convenient for numerical calculations as the LS 
equation. 



II. GENERALIZED QUANTUM DYNAMICS 

Let us briefly review the main features of the formal- 
ism of the GQD developed in Ref.[12]. As is well known, 
the basic concept of the canonical formalism of quantum 
mechanics is that the theory can be formulated in terms 
of vectors of a Hilbert space and operators acting on this 
space. This formalism rests on the postulates, which es- 
tablish the connection between these mathematical ob- 
jects and observables and prescribe how to compute the 
probability of an event. In the canonical formalism they 
are used together with the dynamical postulate accord- 
ing to which the time evolution of a quantum system 
is governed by the Schrodinger equation. At the same 
time, in the Feynman formalism [13] quantum theory is 
formulated in terms of probability amplitudes without re- 
sorting to the vectors and operators acting on a Hilbert 
space. In this approach the following assumption is used 
as the first basic postulate: 

The probability of an event is the absolute square of a 
complex number called the probability amplitude. The 
joint probability amplitude of a time-ordered sequence of 
events is product of the separate probability amplitudes 
of each of these events. The probability amplitude of an 
event which can happen in several different ways is a sum 
of the probability amplitudes for each of these ways. 
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According to this assumption, the probabihty amph- 
tude of an event which can happen in several different 
ways is a sum of contributions from each alternative way. 
In particular, the amplitude < 'il'2\U{t,tQ)\'ipi >, being 
the probability of finding the quantum system in the state 
\^p2 > &t time t, if at time to it was in the state |^i >, can 
be represented as a sum of contributions from all alter- 
native ways of realization of the corresponding evolution 
process. Dividing these alternatives in different classes, 
we can then analyze such a probability amplitude in dif- 
ferent ways. For example, subprocesses with definite in- 
stants of the beginning and end of the interaction in the 
system can be considered as such alternatives. In this 
way the amplitude < ?/'2|C^(t,to)|V'i > can be written in 
the form |l2| 

< ^p2\U{t,to)\^l > = < '(/'2|V'l > 

+ [ dt2 [ dti <i:2\S{t2,h)\^/ji>, (1) 
Jto Jto 

where < '02|'S'(i2, ^i)|'0i > is the probability amplitude 
that if at time ti the system was in the state |'0i >, then 
the interaction in the system will begin at time ti and 
will end at time t2^ and at this time the system will be 
in the state \'ip2 > ■ Here the interaction picture is used. 

According to the above postulate the probability am- 
plitude < V'2|<S'(i2, > can itself be represented by 
the sum of amplitudes for each of the ways in which the 
subprocess with completely specified instants of the be- 
ginning and end of the interaction in a quantum system 
can happen. However, some supplementary assumptions 
about the history of the system are needed. In the Feyn- 
man approach it is assumed that this history can be rep- 
resented by some path in space-time. In this case the 
amplitude < 'ip2\S{t2,ti)\ipi > can be represented by the 
sum of contributions from all paths corresponding to pro- 
cesses in which the interaction begins at ti and ends at ^2- 
If we assume also that the contribution from a single path 
is an exponential whose (imarginary) phase is the clas- 
sical action for this path (the second postulate of Feyn- 
man's theory) and substitute the expression obtained in 
this manner into Eq.(l), we arrive at Feynman's sum- 
over-paths formula for the transitions amplitudes. At 
the same time, in the formalism of the GQD the history 
of a quantum system is represented by the version of the 
time evolution of the system associated with completely 
specified instants of the beginning and end of the inter- 
action in the system. Such a description of the history of 
a system is more general and require no supplementary 
postulates like the second postulate of the Feynman for- 
malism. On the other hand, the probability amplitudes 
< '02|'S'(i2, >, in terms of which the evolution of 
a system is described within the GQD, are used in the 
spirit of Feynman's theory: The probability amplitude of 
any event is represented as a sum of this amplitudes. In 
Ref. [12] it has been shown that the use of the operator 
formalism of the canonical approach allows one to derive 
a relation for the amplitudes < 4'2\S{t2,ti)\ipi > which 
can be regarded as an equation of motion. 



By using the operator formalism, we can represent the 
probability amplitudes < 'ip2\U{t2,ti)\ipi > by the ma- 
trix elements of the evolution operator, which must be 
unitary 

U+it, to)U{t, to) = Uit, to)U+{t, to) = 1, (2) 

and must satisfy the composition law 

U{t, t')U{t', to) = Uit, to), U{to, to) = 1. (3) 

Meanwhile, S{t2,ti) whose matrix elements are < 
ip2\S{t2,ti)\ipi > may be only an operator-valued gen- 
eralized function of ti and t2, since only U{t,to) = 1 + 
/to J^^ dtiS{t2,ti) must be an operator on the Hilbert 

space. Nevertheless, it is convenient to call S{t2,ti) an 
"operator" by using this word in a generalized sense. In 
the case of an isolated system the operator S{t2, ti) can 
be represented in the form 

5(t2,ti) = exp{iHot2)f{t2 - ti)exp{-tHoti). (4) 

As has been shown in Rcf.[|l2|, for the evolution op- 
erator U{t,to) given by (|l]) to be unitary for any times 
to and t, the operator S{t2,ti) must satisfy the following 
equation: 

{t2-ti)Sit2,ti) ^ f dti f dt3 
Jti Jti 

x{ti-t3)S{t2,ti)S{t3,ti). (5) 

A remarkable feature of this equation is that it works as a 
recurrence relation, and allows one to obtain S{t2, ti) for 
any ti and t2, if it is known in an infinitesimal neighbor- 
hood of the point t2 = ti. Since the operators S{t2,ti) 
describe the contributions to the evolution operator from 
the processes in which the interaction in the system be- 
gins at ti and ends at t2, the above means that in order 
to construct the evolution operator it is sufficient to know 
the contributions to this operator from the processes with 
infinitesimal duration time of interaction. It is natural to 
associate these processes with the fundamental interac- 
tion in the system under study. This make it possible 
to use the relation as a dynamical equation. One 
needs only to specify the boundary condition determin- 
ing the behavior of S{t2,ti) in the limit t2 — > ti and 
hence containing the dynamical information about the 
system. Denoting the contribution to the evolution op- 
erator from the processes associated with the fundamen- 
tal interaction by Hint{t2, ti), such a boundary condition 
can be written in the form 

S{t2,ti) ^ H,nt{t2,tl) + o{t'), (6) 

where r = ^2 — ii- The parameter e is determined by 
demanding that Hint(t2,ti) must be so close to the so- 
lution of Eq. (|^) in the limit t2 — > ii that this equation 
has a unique solution having the behavior m) near the 
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point t2 = ti. Within the GQD the operator Hint{t2,ti) 
plays the same role as the interaction Hamiltonian in 
the ordinary formulation of quantum theory: It gener- 
ates the dynamics of a system. Being a generalization of 
the interaction Hamiltonian, this operator is called the 
generalized interaction operator. 

The operator -ffint (^2,^1) describes fundamental pro- 
cesses, starting from which, one can construct the evo- 
lution operator. In this process the system in the state 
lipi > evolves freely up to some time t when, as a result of 
the interaction, the state of the system is jumps abruptly 
into the state |V'2 >, and then the system evolves freely 
again. The contribution from this process into the evo- 
lution operator is of the form S{t2~ti) < ?/'2|>l(ii)|'!/'i >, 
where the delta-function is needed for this contribution 
to be nonzero. Thus in this case the interaction operator 
should be of the form 



H,^t{t2,ti) = S{t2~h)A{ti). 



(7) 



As has been shown in Ref. [12], the dynamical equation 
with the boundary condition given by and is 
equivalent to the Schrodinger equation with the interac- 
tion Hamiltonian Hj{t) = i/2A{t). Thus the dynamics 
governed by Eq. (|^) is equivalent to Hamiltonian dynam- 
ics in the case where the generalized interaction operator 
is of the form 



H,ntit2,ti) ^ -2i6(t2 - ti)Hi{h). 



(8) 



Correspondingly the interaction operator in the 
Schrodinger picture 

Hintit^ - h) = ex.p{-iHot2)Hintit2,ti)exp{iHoti) 



has the form 



(9) 



where Hj — exp{—iHot)Hi{t) cx-p{iHot). In this case 
the interaction generating the dynamics is instantaneous. 
On the other hand, from both the physical and mathe- 
matical points of view there are no reasons to restrict 
ourselves to the case where the interaction operator is 
of the form (^. From the mathematical point of view, 
the boundary condition (||) with the operator Hint{t2, ti) 
given by (^) is not only possible boundary condition for 
Eq.(^. From the physical point of view, this equation is 
a unique consequence of the representation {V) and the 
unitarity condition (|^) . The representation (|l ) in turn is 
a consequence of the first Feynman postulate that, as is 
well known, is formulated as a result of the analysis of 
the phenomenon of the quantum interference and hence 
is one of the most fundamental postulate of quantum 
theory. Thus Eq.(|) is a unique consequence of the first 
principles and can be considered as the most general dy- 
namical equation consistent with the current concepts of 
quantum theory. Note, in this connection, that no new 
fundamental concepts and postulates are used in the for- 
malism of the GQD. A novelty of this formalism consists 



in the fact that some basic postulates of the Feynman 
and canonical approaches to quantum theory are used in 
combination. This allows one to formulate the theory in 
terms of the operators S{t2,ti). As has been shown in 
Ref. [12], theory provides a more detailed description of 
the dynamics of a quantum system then the description 
directly in terms of the evolution operators, as in the case 
of the canonical formalism, or in terms of Feynman's path 
amplitudes. In the case where the interaction operator is 
of the form (^), i.e. the interaction is instantaneous, the 
Schrodinger equation for the evolution operator 



, dU{t, to) 
'' dt 



Hiit)U{t,to) 



and Feynman's sum-over-paths formula follow from the 
representation and Eq.(^). At the same time, the dy- 
namical equation permits the generalization to the case 
where the interaction is nonlocal-in-time, i.e. the time 
durations of the interaction in the fundamental processes 
that determine the dynamics of a system are not zero. In 
this case the dynamics depends not only on the form of 
the operator h\^^{t) but also on its dependence upon 
the duration time t of the interaction. However, as we 
have seen, only the behavior H^^) (r) in the limit r ^ 
is relevant: Knowing the behavior of h['^^{t) in the in- 
finitesimal neighborhood of the point r = is sufficient to 
construct the evolution operator by using Eq.(|). Thus 
within the GQD we deal with a new type of nonlocal- 
ity. In fact, the ordinary way of nonlocalization of a 
quantum field theory consists in introducing a nonlocal 
form factors that depend on parameters determining a 
scale of nonlocality. As for the operator Hint{t2, ti), only 
its values in the infinitesimal neighborhood of the point 
t2 = ti are relevant and hence the scale of its nonlocal- 
ity in time is infinitesimally small. Thus in this case we 
deal with some quasilocal operators. This is very im- 
portant from the point of view of applications to QFT 
where nonlocalization aimed at resolving the problem of 
the UV divergences leads to a loss of covariance of the 
theory. In Ref. [|l^ it has been shown that after renormal- 
ization the dynamics of the three-dimensional theory of 
a neutral scalar field interacting through a (p^ coupling is 
governed by the generalized dynamical equation (^) with 
a nonlocal-in-time interaction operator. 

In order that the dynamical equation (j^) with the 
boundary condition (||) have a unique solution, the oper- 
ator Hi„t(t2,ti) must be sufficiently close to its relevant 
solution. This means that this operator must satisfy the 
condition 



(^2 — ti)Hint(t2, ti) 



X iJ™t {t2 , U)H,nt [h, h) + o(t^+1 ) . (10) 

Note that the value of the parameter e depends on the 
form of the operator Hintit2,ti). Since S{t2,ti) and 
Hint{t2,ti) are only operator- valued distributions, the 
mathematical meaning of the conditions and ( p^ ) 
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needs to be clarified. We will assume that the condition 
means that 

< *2| / dt2 I dhS{h,h)\^i > 



t — to 

21/ dt2 I dtiH„u{t2,h)\^i>+o{T'+^), 



<^2\ I U,i2 

Itn Jt, 







'i'l > and 15*2 > of the Hilbert space. 
10[) has to be considered in the same 



for any vectors 
The condition 
sense. 

If Hint{t2,ti) is specified, Eq.(p[) allows one to find 
the operator ^(i2,<:i). Formula (|l^can then be used to 
construct the evolution operator U {t, to) and accordingly 
the state vector 

\iP{t) \^P{to) > + f dt2 [ dhS{t2,h)\^p{to) >(11) 

J to J to 

at any time t. Thus Eq.(|5|) can be regarded as an equation 
of motion for states of a quantum system. By using (|l|) 
and , the evolution operator can be represented in the 
form 

i f°° 

< n2\U{t,to)\ni >^< n2\ni > +— / dx 

exp[-i(z - En.^)t] exp[i(z - Enjtp] 
{z - £'„J(z - En,) 

X <n2\T{z)\ni>, (12) 

where z = x + iy, y > 0, and 



< n2\T{z)\ni >= i / dr exp{izT) < n2\T{T)\ni > . 
Jq 

(13) 

Here |n > are the eigenvectors of the free Hamiltonian 
Hq, i.e. -ffol^ >= En\n >, and n stands for the entire 
set of discrete and continuous variables characterizing the 
system in full. From (^2|), for the evolution operator in 
the Schrodinger picture, we get 



C/«(t,0) 



2tt 



dxe:iq){-izt)G{z), (14) 



where 



< n2\G{z)\ni >= 



< 7i2|«i > <n2\T{z)\ni> 



z-En, {z - E„^){z - E„J' 

^ (15) 
Eq. (n5[) is the well-known expression establishing the con- 
nection between the evolution operator and the Green 
operator G{z) and can be regarded as the definition of 
the operator G{z). 

The equation of motion (H) is equivalent to the follow- 
ing equation for the T-matrix llj] : 



d < n2\T{z)\ni > 
dz 



E 



< n2\T{z)\n >< n\T{z)\ni > 



{z - En)^ 



with the boundary condition 

<7l2|T(z)|ni> < 7l2|B(z)|71l > +o(|zr'^), 

\z\ — > oo 



where /3 = 1 + e, and 



(17) 



< n2\B{z)\ni >= i / dr exp{izT) < n2\B{T)\ni >, 

B{t) being an arbitrary operator that has the following 
behavior in the limit r — )■ 0: 

< n2\B{r)\n, > - < n2|i?f^)(T)|ni > +o{t^). 

Note in this connection that, while we define the oper- 
ator Hint{t2,ti) for any times ti and t2, only its values 
for infinitesimal duration times t — t2 ^ ti of interac- 
tion are relevant: Knowing the behavior of Hint{t2,ti) 
in the infinitesimal neighborhood of the point t2 — ii 
is sufficient to construct the evolution operator by solv- 
ing Eq.(^. The two interaction operator Hint{t2, ti) and 
Hln,tit2;ti) are dynamically equivalent, provided 



H^nt{t2,tl) - H[^^{t2M)+o{{t2 -hf), 



■0. 



In fact these operators lead to the same solution of 
Eq.(||), i.e. generate the same dynamics. Correspond- 
ingly knowing the behavior of B(z) in the limit \z\ — > oo 
is sufficient for obtaining a unique solution of Eq.(^6|). 
This behavior is uniquely determined by the behavior of 
if|*)(r) in the limit r ^ 0. At the same time, the opera- 
tor H^nti''') ii6ed not be such that the Fourier transform 
c?T exp(izr) < n2\H^^l{T)\ni > exists, since it deter- 
mine the behavior of the operator T{t), for which such a 
Fourier transform must satisfy, only in the limit r — > 0. 
Of course, one can use a dynamical equivalent operator 
for which such a Fourier transform exists. The opera- 
tors B(t) is an example of such operators. However, it 
is not convenient to deal with such interaction operators. 
In fact, in this case one has to take care of the behav- 
ior of ff|^)(r) not only in the limit t ^ but also in 
the limit r — > oo. Nevertheless, formally we can con- 
struct the operator B{z) for any z, by using the operator 
B{t) that being dynamically equivalent to the operator 

hI^^{t) satisfies the above requirement. 

It should be noted that the T-matrix obtained by solv- 
ing Eq.(|l^) satisfies the following equation: 



< n2\T{zi)\ni > - < n2\T{z2)\ni >- 



{z2 - zi) 



^ < n2\T{z2)\n^ >< n\T{zi)\ni > ^^^^ 



{Z2 - En){zi - En) 



(16) 



Below we will show that in some cases it is convenient 
to use this equation itself for obtaining the T-matrix. 
Note also that in general the interaction operator has 
the following form [ p^ : 

H,nt{t2M) = - tl)Hi{ti) ^ Hnon{t2M), 
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where the first term on the right-hand side of this equa- 
tion describes the instantaneous component of the inter- 
action generating the dynamics of a quantum system, 
while the term Hnonit2,ti) represents its nonlocal- in- 
time component. 

III. NUCLEON DYNAMICS IN LEADING 
ORDER OF THE EFT APPROACH 

Let us now show that the leading order nucleon dynam- 
ics in the effective theory of nuclear forces is governed by 
the generalized dynamical equation (||) with a nonlocal- 
in-time interaction operator. For the sake of simplicity 
we will consider the nucleons as identical, spinless par- 
ticles (describe by a field ■0) with 3-momenta Q much 
smaller than their mass m, the mass difference A to their 



first excited state, the pion mass /i, and the range of 
their interaction. Being straightforward, generalization 
to spin and nonidentical particles cannot give rise to an 
essential change of the dynamical situation in a nucleon 
system. Thus we will consider the dynamics of spinless 
particles at wave lengths which are large compared to 
the range of their interaction. Such a dynamics can be 
described by the EFT of short-range forces developed by 
van Kolck [Q. In this theory particles are described by 
a field ip, and the effective Lagrangian involve arbitrary 
complicated operators of only ip and ifs derivatives. It is 
assumed these derivatives are associated with factors of 
I/to, 1/A, or and, therefore, that the effective La- 
grangian can be written as an expansion in c}/(to. A, fi). 

The effective Lagrangian of the theory can be written 
as |l| 



2m 8to-^ 2 8 

_^)ijf . ^r + (^ - ^)ijf - * + ** + (^ - ^f-^} + -{C2- C01' + *V2(5' + *) + . . . , (19) 



where the C2„s are parameters that depend on the details 
of the dynamics of range ~ 1/M, where M characterizes 
the typical scale of all higher-energy effects. Canonical 
quantization leads to familiar Feynman rules, and the ^ 
propagator at four-momentum p is given by 



^(/,P) 



„o _ Pi , _Pl. 

2m ' 8m3 



- le 



(20) 



The four- 5" contact interaction is given by —iv{p,p'), 
with 

v{p,p') = Co + C2(p2 + p'') + 2C^p •?' + ..., (21) 

p(p') being the relative momentum of the incoming (out- 
going) particles. The problem can also be easily solved 
by using the time-ordered perturbation theory, since in 
this case we deal only with the particles evolving forward 
in time. Obviously this way is more convenient for con- 
structing the off-shell T-matrix. 

Let us consider the two-particle system at energy E = 
.2 1 4 

- — 7*^- + ... in the center-of-mass frame. The key point 
of the EFT approach is that the problem can be solved 
by expanding in the number of derivatives at the vertices 
or particle lines. In leading order one has to keep only 
the first term Co in (^^7 and correspondingly only the 
first two terms p'^ and ^ in the propagator ( pO| ) . In this 
case the two particles evolve according to the familiar 
nonrelativistic Schrodinger propagator 



Go{z) = 



d^k |k >< k| 
(27r)3 z _ + J, 



(22) 



In this order the Lagrangian can be rewritten in the form 

£ = 1'+(i9o + — V^)* - ico^'+**+^'. (23) 
2to 2 

Conservation of particle number reduces the two-nucleon 
T-matrix to a sum of bubble diagrams. The ultravio- 
let divergences can all be absorbed in the renormalized 

(ft) 

parameter Cq . Summing the bubbles to a geometric 
series, one gets the T-matrix p4] 



<P2|r(o)(z)|pi >=- 



47r 



(24) 



What is interesting the T-matrix given by ( p4| ) is ex- 
actly the same as in the model [12,16]. This model pro- 
vides an example showing that the GQD allows one to 
extend quantum dynamics to the case where the interac- 
tion generating the dynamics of quantum system is non- 
local in time. Thus the model shows that the situation 
where the dynamics of a quantum system is generated 
by nonlocal-in-time interaction is possible in principle. 
From the above it follows that this possibility is realized 
in low energy nucleon dynamics. The same T-matrix can 
be obtained via the LS equation with the potential 



T/(p2,Pi) = 



(25) 



by using some regularization and renormalization proce- 
dures. In order to show this fact, let us consider the 
evolution problem for two nonrelativistic particles in the 
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c.m.s. Assume that the generalized interaction operator 
in the Schrodinger picture is of the form 



<P2|i?i^l(T)|pi >=^*(p2Mpi)/(r), 



(26) 



where /(r) is some function of the duration time r of 
interaction, and the form factor (^(p) must have the fol- 
lowing asymptotic behavior for |p| — > oo: 



'P{p)-\pr, (IpI-oo). 



(27) 



If the interaction is instantaneous, then /(r) — — 2iA(5(T), 
and 

< P2|ffit)(T)|Pi >= -2iS(h - h) < P2|y|pi >, (28) 
< P2|^^|pi > being the separable potential 

<P2|l/|pi >= A^*(p2)¥'(pi). (29) 
The solution of the LS equation with this potential is 

<P2|r(z)|pi >=¥>*(p2)<^(Pl) 
d^k |(^(k)|2 



(2^)3 {z - Ek) 



(30) 



In the case a < ^, the potential (25) does not make 
sense without renormalization, since it gives rise to UV 
divergences. In fact, in this case the integral in (30) is 
not convergent, and one has to use some regularization 
procedure. Using the dimensional regularization we can 
write the following expression for the regulated T-matrix: 



< P2|rg(z)|pi (^*(p2)</3(Pl) 



(2^)3 {z-Ek), 
This expression can be rewritten in the form 

<P2|Te(z)|pi >=(p*(p2)^(pi) 

1 f dp-^k |<^(k)|2 



(31) 



where (C^^))-! = A^^ - / 



(27r)3 (z - Ek)Ek 



(32) 



(27r)3 



Let us use the 



scattering length ao as an additional experimental infor- 
mation that is needed for renormalization. The parame- 
ter Cg'^'' in the expression (32) for the T-matrix is related 
to the scattering length by 



The above means that we have to fix the value c'^'^ which 
we will denote by Cq^\ Now we can let e — ^ in Eq.(32), 
and get 



< P2|T(2:)|pi >= (/3*(p2)<^(Pl) 



dPk |^(k)|2 



(R) 



(2^)3 (z - Ek)Ek 



(33) 



It is easy to see that this renormalized T-matrix is not 
a solution of the LS equation with some potential be- 
cause of the slow rate of decay of the form factor fip) as 
p — > oo. Note in this connection that the renormalized 

.(e) r d'^-'k |y(k)lS-i 



strength of the potential Ae = {Cq' — J -(2tt)^ 



-r 



tends to zero as e — s- 0. This means that after renor- 
malization the interaction Hamiltonian is zero. This is 
a manifestation of the well-known fact that one cannot 
construct a renormalized Hamiltonian, and there are not 
any equations for renormalized amplitudes within Hamil- 
tonian formalism. At the same time, one can easily verify 
that the T-matrix given by (33) satisfies the generalized 
dynamical equation (||). Moreover, as stated above, it 
coincides with the T-matrix of the model [12,16]. In this 
model it is assumed that the generalized interaction oper- 
ator is of the form (26). In this case, < p2\S{t2,ti)\pi > 
can be represented in the form 

< P2|5(t2,tl)|pi >- ^*{p2MPl)S{t2,tl). (34) 

Correspondingly, the T-matrix should be of the form 

< P2|T(z)|pi >= (p*(p2)(^(pi)t(z), (35) 

From (16) it follows that the function t{z) must satisfies 
the equation 



dt{z) ^ f d^k |^(k)|2 



dz " J (27r)3 (z - Ek)^ ' 

This equation can be rewritten in the form 

di{z) f d^k |(p(k)|2 



dz 



(2^)3 (z - i?fc)2 



(36) 



(37) 



with t(z) = (i(z)) . All solutions of Eq.(^6|) satisfying 
the condition 



t*{z) ^ t{z), z e (-00,0), 



(38) 



are physically realizable [12]. Each of them corresponds 
to the definite function /(r) in the generalized interaction 
operator (26). As we have noted, only the asymptotic be- 
havior of the function /(r) as r — > is relevant, and for- 
mally one can use any function that has this asymptotic 
behavior. Knowing the asymptotic behavior of H^^^^ir) 
in the limit r — > allows one to determine the asymp- 
totic behavior of the function t{z) in the limit z —^ oo 
(the one-to-one correspondence between these behaviors 
follows from the Fourier transform). This behavior in 
turn can be used as a boundary condition for the differ- 
ential equation (36). The T-matrix obtained in this way 
can then be used for constructing the evolution operator. 
Thus knowing the asymptotic behavior of h\^^(t) in the 
limit T ^ uniquely determines the dynamics of the sys- 
tem. On the other hand, as it follows from Eq.(37), the 
asymptotic behavior of t{z) is determined by the large- 
momentum behavior of the form factor </'(p). In the case 
a > 5, at it is easily seen from Eq.(37), all solutions of 



8 



Eq.(^6|) tend to nonzero constants as 
have 



tiz) 



A. 



(39) 



This means that in this case the only possible form of the 
function /(r) is 



/(T) = -2*A5(r)+/'(r), 



(40) 



where the function /'(t) has no such a singularity at the 
point r = as the delta-function. Thus, in the case 



a > 



the interaction operator should be of the form 



< P2|i?i:l(r)|pi >= -2td{T) < p2|^|pi >, (41) 

< P2|V^|Pi > being the separable potential 

<P2|^^|pi A^*(p2)^(pi) (42) 

and hence the dynamics generated by this operator is 
equivalent to the dynamics governed by the Schrodinger 
equation with the separable potential Ay?* (p2)v?(pi). 
Solving Eq.(|3^) with the boundary condition (39), for 
the T- matrix, we easily get the expression (30), i.e. in 
the case a > i, the interaction in the system may be 
only instantaneous. 

Let us now consider the case — | < a < | (the re- 
striction a > — ^ is necessary for the integral in (^) 
to be convergent) where the ordinary separable-potential 
model does not make a sense without renormalization. In 
this case solutions of Eq.(|3q) tends to zero as \z\ ^ oo: 



where bi = —4TTcos{an)m'^~i, and 62 is some arbitrary 
constant. Thus the T-matrix has the following asymp- 
totic behavior for Izl — > cxi: 



<P2|T(z)|pi > 

(bi{~zr-i +b2{-z) 



2q- 



V*(P2)V(P1 



|2q-1- 



(43) 



By using asymptotic methods (see, for example, Ref.Q) 
it is easy to show that this large z behavior of the function 
t{z) corresponds to the following behavior of T(t) related 
to the T-matrix by Eq.(13) in the limit t — > 0: 



<P2|T(r)|pi> ^ifi* {P2MP1) 



X aiT 



a2T 



+0(t-2") 



(44) 



where ai = —ibiT ^{^ — a)exp[i{—^ + j)n], and 
0-2 = ^2r^^(l — 2a)exp{—iaTT), T{z) being the gamma- 
function. On the other hand, according to Eqs.(4) and 
(6) the asymptotic behavior of < P2|T(t)|pi > is deter- 
mined by the interaction operator -fff*) (t) 

< P2|r(r)|pi > - < P2|4:)(t)|pi > +o(r^), (45) 



00, i.e. we where the parameter e is determined by demanding that 
exp(i_ffoi2)^fi,^](T) exp(— iifoti) must be so close to the 
relevant solution of Eq.(|^) in the limit r — > that this 
solution is unique having the asymptotic behavior (45). 
It is easy to see that in the separable case all solutions 
of Eq.(g6|) for given a have the same leading term in the 
expansion (43), since the parameter bi is uniquely deter- 
mined by the value of the parameter a. Only the second 
term containing free parameter 62 distinguishes the dif- 
ferent solutions of this equation. On the other hand, the 
first two term in the expansion (44 ) are uniquely deter- 
mined by the first two term of the expansion (43). This 
means that the parameter e in (45) is equal to 2a, and 
hence the generalized interaction operator H^^] (r) should 
be of the form 



<P'l^)l 



r)|p >= ip 



a2T 



(46) 



Of course, the interaction operator (^6|) may be supple- 
mented by any term being of order o(t~^"). However, 
this will not change the solution of Eq.(||), i.e. in the 
is case we will get the generalized interaction operator 
being dynamically equivalent to the interaction operator 
(^). As we have noted, the operator H^^^{t) describes 
the contributions to the evolution operator from the pro- 
cesses with infinitesimal duration time r of interaction, 
and only its behavior in the limit r — s- is relevant for 
us. 

The solution of Eq.(37) with the initial condition t{z ~ 
0) = to is 



i{z) ^ia- z 



(27r)3 {z - Ek)Ek 



(47) 



From this we easily get the expression (33) for the T- 
matrix. This solution of Eq.(16) has the behavior (43), 
provided — ~^2^r^- Taking this fact into account, one 
can rewrite this solution in terms of the parameters hi 
and 62 



< P2|r(z)|pi >= Niz)ip*{p2Mpi), (48) 



where 



(49) 



with 



f |y(k)p-|k|-^° 

By using Eqs.(12) and (47), we can construct the evo- 
lution operator 



(51) 



< p2|C/(t,io)|Pi >=< P2IP1 > +7^ / dx 

J-00 

eyip{-izt)N{z)ip*(p2)if{pi) 
{z - Ep^){z - Ep,) 
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where z = x + iy, y > 0. The evolution operator (|51^ 
satisfies the composition law and is unitary provided 
the parameter 62 is real. Thus, in the case — ^ < a < |, 
the dynamics is well defined, and all the requirements of 
quantum theory are satisfied. 

The model under consideration is an transparent il- 
lustration of the possibility of the extension of quantum 
dynamics provided by the GQD. From the point of view 
of Eq.(36), the restriction to Hamiltonian dynamics is 
equivalent to the restriction to the solutions that tend to 
some nonzero constant in the limit |z| — )■ 00. In this case 
the interaction in the system is instantaneous, and the 
dynamics is Hamiltonian. At the same time, this equa- 
tion has solutions that tend to zero as |2:|— s-oo. This 
takes place in the case — ^ < a < 5. In this case the 
interaction is nonlocal in time and is parameterized by 
the interaction operator (^). In this case the dynamics 
is non-Hamiltonian. There are no reasons to restrict our- 
selves to the solutions of Eq. (^) which tend to nonzero 
constants, because all solutions of this equation satisfy- 
ing the condition (38) are physically realizable, and such 
an extension of quantum dynamics is possible in princi- 
ple. The remarkable fact is that this possibility is realized 
in low energy nucleon dynamics, and, as we have noted, 
the leading order contribution to the nucleon-nucleon T- 
matrix is described by our model. In fact, in the case 
■^(p) — 1, for the T-matrix given by (47), we have 



<P2|T(z)|pi >= 



62 - ^bl^/z 



3 

11712 , 



^0 



(R) 



where bi 



m\/rn 



, and Cq 



{R) 



b2 



(52) 



In this case the 



generalized interaction operator is of the form 



<P2|i?i:)(r)|pi >= 



4-y/7r exp(i 



167r2 



(53) 



This operator parametrizes the leading order contact 
component of the NN interaction and generates the two- 
nucleon dynamics that is described by the evolution op- 
erator 



< p2|C/(t,io)|Pi >=< P2IP1 > 



27r 



dx 



exp{—izt)N{z) 
(z - Ep^ ){z - Ep, 



(54) 



It should be noted that the interaction operator (|5S 
tains the renormalize parameter Cq^\ while in the or- 
dinary methods based on the use of a renormalization 
procedure this parameter appears in renormalized am- 
plitudes at the final stage. The generalized dynamical 
equation (^ with the interaction operator (53) is well 
defined and allows one to describe, in a consistent way, 
the dynamics generated by the leading order component 
of the NN interaction. 



IV. NON-HAMILTONIAN CHARACTER OF 
LOW ENERGY NUCLEON DYNAMICS 

Let us now discuss advantages of the GQD and new 
possibilities that it opens for describing low energy nu- 
cleon dynamics. As we have seen, the T-matrix (33) ob- 
tained by solving the dynamical equation (16) with the 
interaction operator (46), can be also obtained starting 
with the singular potential 

l/(p',p) = A^*(p')<^(p), <^(p)'-|pr", a<i. 

However, in this way one cannot determine a potential 
that could parametrize the interaction in the system. The 
above singular potential does not make sense without 
renormalization and is only of formal importance for the 
problem under consideration. In fact, as we have shown, 
the strength of potential, which in the regulated LS equa- 
tion is adjusted to give the correct scattering length, be- 
comes zero after removing regularization. All the infor- 
mation contained in this potential is that the T-matrix 
is of the form (35) with the same form factor. This in- 
formation is not sufficient to construct the T-matrix. In 
addition, one needs to use some experimental data, for 
example, the scattering lengths. However, one cannot 
construct an interaction Hamiltonian that could contain 
all this information. The T-matrix (33) obtained in this 
way does not satisfy the LS equation. Moreover, as we 
will show below, it has the properties that are at variance 
with the Hamiltonian formalism. Thus in this case we 
have only a calculation rule that allows one to compute 
the T-matrix: One cannot derive any renormalized equa- 
tion containing a well-defined potential that parametrizes 
the interaction in the system. This problem is the price 
for trying to describe the dynamics of the system af- 
ter renormalization in terms of the Hamiltonian formal- 
ism, despite this dynamics is non-Hamiltonian. From the 
more general point of view provided by the GQD we see 
that the T-matrix (33) satisfies the generalized dynamical 
equation (16) with the nonlocal-in-time interaction oper- 
ator (46), and this operator describes the fundamental 
interaction in the system. Knowing this operator is suffi- 
cient to construct the T-matrix and hence the evolution 
operator, and in this case we deal with the well-defined 
equation that does not require renormalization. As we 
have seen, an example of such a dynamics is the leading 
order nucleon dynamics where, as a consequence of the 
symetries of the underlying theory, the form factor must 
satisfy the condition (p{p) — 1. 

In order to clarify the character of the leading order 
nucleon dynamics let us examine the properties of the 
evolution operator of the more general theory with the 
interaction operator (46). In the Schrodinger picture, 
this operator < p2|y(i)|pi >=< p2|C/s(t, 0)|pi > can be 
rewritten in the form 



< P2|V'(t)|pi >=< P2IP1 > e^p{~iEp^t) 
i exp(-izt) <p2|r(z)|pi > 

SttI^ {z - Ep,)iz - EpJ ' 



(55) 
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where < p2|T(z)|pi > is given by (^8|). Since this T- 
matrix satisfies Eqs.(^6|) and (|l8|), the evolution oper- 
ator (51) is unitary, and satisfies the composition law 
(^. Correspondingly the operators V{t) constitute a one- 
parameter group of unitary operators, with the group 
property 



V{h+t2)=V{h)V{t2), V{0) 



(56) 



Assume that this group has a self-adjoint infinitesimal 
generator H which in the Hamiltonian formalism is iden- 
tified with the total Hamiltonian. Then for >S 'D{H) 
we have 



V{t)\ij>^\ij> ^ _ 
t t^o 

From this and (55) it follows that 



(57) 



The above means that low energy nucleon dynamics 
in the effective theory of nuclear forces cannot be gov- 
erned by the Schrodinger equation and hence is non- 
Hamiltonian. Correspondingly the interaction of nucle- 
ons cannot be parametrized by an interaction Hamilto- 
nian defined on the Hilbert space. At the same time, as 
has been shown on the example of the effective theory 
with the Lagrangian (23), this dynamics is governed by 
the generalized dynamical equation (H) with a nonlocal- 
in-time interaction operator. In the above theory describ- 
ing the leading order nucleon dynamics this operator is of 
the form (53), and two nucleon dynamics is described by 
the evolution operator (54). This evolution operator is 
unitary and satisfies the composition law (|^). Thus the 
use of Eq.(|^) as an equation of motion in the effective 
theory of nuclear forces allows one not only to calculate 
scattering amplitudes but also to construct the evolution 
operator describing low energy nucleon dynamics. 



with 



< P2|i?/|pi >= 



I 

2^ 



^ .<P,|r(.)|p > 



where z — x + iy, and y > 0. Since < p2|r(2;)|pi > 
is an analytic function of z and, in the case a < i, 
tends to zero as |z|^oo, from Eq.(58) it follows that 
< P2|^^/|Pi for any p2 and pi, and hence H — Hq. 
This means that, if the infinitesimal generator of the 
group of the operators V{t) exists, then it coincides with 
the free Hamiltonian, and the evolution operator is of the 
form V(t) — exp(— liJoO- Thus, since this, obviously, is 
not true, the group of the operators V{t) has no infinites- 
imal generator, and hence the dynamics is not governed 
by the Schrodinger equation. 

It should be also noted that in the case a < i, S{t2, ti) 
is not an operator on the Hilbert space. In fact, the wave 
function 



V'(p) =< pIV- >=< p\s{t2,h)\^Ji >-- 
(fik 



(/7*(p)l(i2,<i 



(59) 



is not square integrable for any nonzero I'i/'i >j because of 
the slow rate of decay of the form factor y'(p) as |p| oo. 
Correspondingly, in the case a < the T-matrix given 
by (48) does not represent an operator on the Hilbert 
space. However, as we have stated, in general S{t2,ti) 
may be only an operator- valued generalized function such 
that the evolution operator is an operator on the Hilbert 
space. Correspondingly the T-matrix need not be an 
operator on the Hilbert space. It is enough that the evo- 
lution operator given by (51) is such an operator. The 
T-matrix and 5(^2, ^i) satisfy these requirements not only 
for a > ^ but also for — ^ < a < ^, since the evolution 
operator (51) is an operator on the Hilbert space. At the 
same time, in the case a < ^ we go beyond Hamiltonian 
dynamics. 



V. NEXT-TO-LEADING ORDER 

We have shown that in leading order the dynamics of 
the EFT under consideration is governed by the gener- 
alized dynamical equation with a nonlocal-in-time inter- 
action operator. Let us now show that this is true also 
in next-to-leading order. For this let us come back to 
the theory with the Lagrangian (19) in which it is as- 
sumed that particles evolve forward in time, and particle 
number is conserved. Summing the graphs of the pertur- 
bation series yield the following on-shell T-matrix fl^ : 



T„,(fc,p'.p) = T(°)(fc) 



2C2^^^A:p' • p 



(60) 



where Tis"* (fc) is the S-wave amplitude and is given by 



C. 



(R) 



link 



G 



{R.) 



(cry 



1 



fc2 



x[l + 0((Q/Af)4)] 



(61) 



^ = |Pi| = IP2I, Pi/k, Cq^\ €2^^ and Cj^"'' are renor- 
malized parameters. In this way one can also con- 
struct the off-shell T-matrix. At the same time, as it 
has been shown in leading order, the T-matrix can be 
obtained without summing the diagrams and resorting 
to regularization and renormalization procedures, pro- 
vided that the UV behavior of the T-matrix elements 
as functions of momenta is known. As we have seen, 
the requirement that the T-matrix being of the form 
< P2|r(z)|pi >= t{z) satisfies the generalized dynam- 
ical equation (16) yields the formula (52) where only the 
parameter Cq is free. In next-to-leading order the re- 
lation between the Lagrangian (19) and the form of the 
T-matrix is not so straightforward, because in this case 
one has to take into account the fact that this Lagrangian 
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is only of the formal importance, and should be supple- I 0[ (QA^)^] (62) 

mented by some counterterms. Nevertheless, let us as- mM 

sume that they do not effect on the form of the T-matrix. 

In this case from the analysis of the perturbation series 

of the theory it follows that the two-particle T-matrix (in 

this paper we focus on the S-wave channel) should be of 

the form 

where Ci = Q^VQ^^ Substituting (|6|) into (16), we 
<P2|T(z)|pi >= ti{z) + CM^){pI + P2^) + get 



dti{z) 
dz 



d^k 1 



t?(z) 1 + 



4m3(z - !^) 
I 



+ 0[^iQ/M)% (63) 



Since the problem with UV divergences must not arise in Correspondingly, the function i{z) = {t{z)) satisfies 
treating the generalized dynamical equation, we have the following equation: 



1 4j^2 ^ 1 



(64) 



with 6Ci = 0{j^) and 



t,(z)^t,iz)^t{z)^l + 0[i^r]], (65) 
where t{z) is a solution of the equation 



di{z) f d^k ( 2zfc2 

J (2^)3(z-^)2 l^+2m2(z-^) 



(67) 



dt{z) 
dz 



d^k 



2zP 



(2^)3(z-^)2 \ ■ 2m^{z-^) 



fe2 - 



(66) Its solution with the initial condition i{z = 0) — to is 



i{z) = io + / ds 



d^k 



2sk^ 



(27r)3 (s- fei)2 I 2m2(s- ■^) 



to 



d'k z+^ z^ 

(2-)^ 4m(--^)^(-^). 



zr7i2 ^ / z \ 
to-—V^{l + I 



An 



8m J 



Denoting = thj, for the function t{z), we get 



(68) 



t{z) 



1 



11712 



MR) ■ Att V ' 8m. 
Substituting (65) with t{z) given by (69) into (63) yields 



V 8m/ 



[1 + 0{{Q/M)% 



<P2|r(z)|pi >= — J 



1 + ^(p?+p1) 4. 



(69) 



(70) 



One can obtain the on-shell T-matrix, by putting |pi| = |p2| — k, and z = ^ — 



Tosik,p'-p) = - 1 + 2 



C 



^0 



^0 



fc4 



(R) 4tt \ m Arrfi 



8m2 32m4 
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+0[-±j(Q/M).| 



a 



(R) 



An 



2^ J 



(71) 



Thus the requirement that the T-matrix of the form ( p^ ) 
satisfies the generahzed dynamical equation (16) yields 
the formula (70), and in this way we get the same two- 
particle on-shell T-matrix that in Ref.|Q has been ob- 
tained by summing the bubble graphs of the EFT. 

In order to obtain the form of the generalized interac- 
tion operator that leads to the T-matrix (70), one has to 
examine the large z behavior of the function t{z) in the 
limit \z\ oo. From (69) it follows that 



hi 



o[\s\ 



(72) 



where s = 



(1+4^) 



hi 



47r 



and hf 



Of course, the asymptotic behavior of t{z) can be repre- 
sented in the form of the expansion in z^^: 



t{z) 



b[ 8mb[ 



(-Z)3 (-z)t i-z)' 



+ oi\zr% (73) 



_ (fell! 

- cr 



However, for the prob- 



where K = ^ and foo 

lem under consideration the values of z much larger than 
Q but much smaller than m should be considered as in- 
finitely large, and the behavior of the function in this 
region is described by the expansion (|7^). The interac- 
tion operator H^^I{t) that gives rise to this behavior of 
t{z) is 



<P2|i?i:l(r)|pi >= 1 



C. 



(R) 



(R) 



4V7Fexp(i^) 167r2 



{R) 



(74) 



8m 



This operator parametrizes the contact term of the NN 
interaction in the S-wave channel up to next-to-leading 
order. The generalized dynamical equation with this in- 
teraction operator uniquely determines the dynamics of 
the system. The equation of motion (|^) with this inter- 
action operator is well-defined and allows one to obtain 
the T-matrix and the evolution operator without resort- 
ing to regularization and rcnormalization. It is easy to 
show that the evolution operator constructed in this way 
(the evolution operator (12) with the T-matrix given by 
(70)) is unitary and satisfies the composition law (||) up 
to next-to-leading order. 



VI. INTEGRAL EQUATIONS 

As we have shown up to next-to-leading order of the 
EFT with the Lagrangian (19), the dynamics of nucleons 



at low energies is governed by the generalized dynamical 
equation (^ with the nonlocal-in-time interaction oper- 
ator (fz^), provided the parameters of the theory satisfy 
the condition (64). Within the GQD this operator is well 
defined. Of course, the dynamical situation in the theory 
with such an interaction differs from that in the theory 
with ordinary potentials. The operator5'(t2, ii) is only 
an operator-valued generalized function on the Hilbert 
space. However, as we have noted this is not at variance 
with the general requirements of quantum theory, since 
the evolution operator given by (51) is an operator on 
this space. The generalized dynamical equation (^) is not 
equivalent to the Schrodinger equation in this case. Nev- 
ertheless, it is a well defined equation, and allows one to 
obtain the T-matrix without resorting to regularization 
and renormalization procedures. This is very important 
for practical calculations, since for solving realistic prob- 
lems one has to deal not only with the contact component 
of the NN interaction but also with its long-range one. 
In this case the interaction operator is of the form 



< p2|i?i:l(T)|pi >= -2*<5(r)F(p2,Pi) 

+ < P2|^^non(T)|Pl >, 



(75) 



where < P2\Hnon{T)\pi > represents the contact 
nonlocal-in-time component and V^(p2,Pi) is a potential 
describing the long-range component of the NN interac- 
tion. In general it consists of the meson-exchange po- 
tentials and the Coulomb potential in the proton-proton 
channel. In the leading order the nonlocal component is 
given by (53) and for the interaction operator, we can 
write 



<P2|i/i:)(r)|pi >= 



4V^exp(i^) 167r2 



7712 VT 



73(7, 



(R) 



-2zd{T)V{p2,Pi). (76) 



In the case of such an interaction operator, the solution 
of the dynamical equation (16) can be represented (see 
Appendix A) in the form 

< p2|T(z)|pi to{z) + ti{z; pi) + ti{z; P2) 

-)-t2(z;p2,Pi), (77) 

where t2{z;p2,Pi) is a solution of the equation 
i2(^;;P2,Pi) = V"(p2,pi) 



d^q K{z;p2,q) , . 
+ I T:r-r^ — <2(2;q,Pi), 



(27r)3 



(78) 



with 



, . f d^k F(p2,k) , , , , 
if(z;P2,q)=A^(z) j (2^)3 ^ \ +^(P2:q): (79) 
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and the functions to(^) E^nd p) are defined as 

(fk t2(z;p, q) 



ii(z;p) ==iV(z) 



io(z) = N{z) 1 + 



(2^)3 z-Ek ' 



d^g ti(z;q) 
(27r)3 z-Eq 



(80) 



where 



1 im^ yfz 



MR) At: 

^^0 



Equations with the next-to-leading order corrections can 
be derived in the same way. 

In Weinberg's power counting the one-pion-exchange 
potential is of leading order. Hence in this order the NN 
interaction operator can be expressed as 



<P2|4:i(r)|pi > = 



4Y^exp( 



;37r^ 



m2 ^ 



(R) 



-2zJ(r)F,(p2,pi), (82) 



where V^(p2,Pi) is the conventional one-pion-exchange 
potential. Substituting this potential into Eq.(|82|) and 
solving it numerically, one can easily obtain the T- matrix 
and hence the evolution operator. Note that conventional 
way of solving the above problem is the formal use of the 
potential 



Vb(p2,pi) = c + y^(p2,pi) 



(83) 



(see, for example, Refs.[7,ll]). We say "formal" since the 
use of such a potential leads to UV divergences, and the 
Schrodinger and LS equations require regularization and 
renormalization. On the other hand, as we have shown, 
the contact interaction, which in (83) is formally repre- 
sented by the term C, is parametrized by the operator 
(53) (the first two terms in the operator (82)). In this 
case we deal with the well defined interaction operators 
and Eq.(^) which does not require regularization and 
renormalization. By using Eq.(57), one can obtain the 
T-matrix so easily as in the case of the pure one-pion- 
exchange potential. 

In order to take into account electromagnetic correc- 
tions, one has to include the Coulomb potential into the 
NN interaction operator 



<P2ii/f:i(T)ipi >= 



4y¥exp(if ) 



m2 ^ 



16^2 



(R) 



-2z<5(r)K(P2,Pi)-2z,5(T 



ei2Q! 



(84) 



where ei2 = 1 for the proton-proton channel and zero 
otherwise. Eq.(75) with the interaction operator (Q) can 
be easily solved numerically. Thus this equation allows 
one to investigate electromagnetic corrections without re- 
sorting to renormalization. 



VII. SUMMARY AND DISCUSSION 

We have shown that the formalism of the GQD al- 
lows one to formulate the effective field theory of nuclear 
forces as an internally consistent theory with equations 
that do not require regularization and renormalization. 
It has been shown that the effective NN interaction is 
nonlocal in time, and low energy nucleon dynamics is 
governed by the generalized dynamical equation (^) with 
a nonlocal-in-time interaction operator. In leading order 
of the EFT approach this operator is given by (p3|). It 
should be noted that the generalized interaction opera- 
tor Hint{t2i ii) is not an effective interaction Hamiltonian 
(only in the local case they are related by Eq.(8)), and, in 
contrast with such Hamiltonians, the use of the general- 
ized interaction operators permits a natural parametriza- 
tion of the NN interaction: The generalized interaction 
operator ( |53| ) generates the unitary evolution of a nucleon 
system. In Sec. Ill we have demonstrated the advantages 
of the GQD in describing the leading order nucleon dy- 
namics in comparison with the ordinary methods based 
on the use of the singular potential < P2|V^|pi >= Cq. 
This potential does not make sense without renormal- 
ization and hence does not contain all the needed dy- 
namical information. In addition, one needs to use some 
empirical data as renormalization constants. In contrast 
with this singular potential, the generalized interaction 
operator (53) contains all the needed dynamical informa- 
tion. For example, it contains the observable, renormal- 
ized parameter Cq . The dynamical equation with the 
interaction operator (^3|) is well defined and allows one to 
construct the T-matrix and the evolution operator with- 
out resorting to regularization and renormalization pro- 
cedures. The dynamical information contains not only in 
the form of the operator ( p3| ) that, as has been shown, is 
a consequence of the symmetries of QCD but also in its 
dependence on the time duration r of the interaction in 
a system. 

As has been shown, only values of H^^^I{t) in the in- 
finitesimal neighborhood of the point r = 0, i.e. at scales 
of the underlying theory are relevant. As already stated, 
the operator H^^^{t) describes the contribution to the 
evolution operator from the processes with infinitesimal 
time duration of interaction. The above means that 



<^2|4:)(t)iv^i > 



<^2|5,(t2,il)|V'l > (85) 



where < '02|'5g(i2, ii)|V'i > describe the contributions to 
the evolution operator from the processes in which the 
quark and gluon degrees of freedom can come to play. 
Here one of the advantages of the formulation of the the- 
ory in terms of the operators S(t2, ti) becomes apparent: 
The relevant amplitudes < ^p2\S{t2,tl)\^pl > of the un- 
derlying theory can be directly used as the matrix ele- 
ments of the generalized interaction operator Hint{t2,ti) 
generating low energy dynamics. These degrees of free- 
dom manifest themselves through the r dependence of 
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the operator (|5^). Thus the parameter in Eq.(|5^) 

parametrizes the leading order effects of the underlying 
physics on low energy dynamics, and hence can be com- 
puted in terms of parameters in QCD. At the same time, 
this parameter can be determined from low energy ex- 
periments. 

The evolution operator governing the two-nucleon dy- 
namics in leading order is of the form (54). In Sec. IV 
we have shown that the group of these operators has no 
infinitesimal generator and hence the dynamics is non- 
Hamiltonian. This means that within the EFT approach 
there are no potential satisfying the requirements of or- 
dinary quantum mechanics that could govern low energy 
nuclcon dynamics. For example, the T-matrix given by 
(52) does not satisfy the LS equation and has properties 
that are at variance with the Hamiltonian formalism. On 
the other hand, as we have seen, the T-matrix satisfies 
the generalized dynamical equation (16) with the interac- 
tion operator given by (|3|). The essential lesson we have 
learned from the above analysis is that many problems 
of the EFT of nuclear forces arise because of ignoring the 
fact that low energy dynamics produced by this theory 
is non-Hamiltonian. For example, the above statement 
that there are not any equations for renormalized ampli- 
tudes in a subtractive EFT should mean that there are 
no such equations within Hamiltonian formalism. From 
a more general point of view provided by the GQD we see 
that such equations exist, and the theory can be formu- 
lated in a completely consistent way. An EFT describes 
low energy physics in terms of a few parameters, and 
these low energy parameters can be computed in terms 
of a more fundamental high energy theory. It is remark- 
able that the interaction operator (|5^) parametrizing the 
leading order contact component of the NN interaction is 
uniquely determined by the symmetries of QCD and con- 
tains one of these parameters Cg^' . Other renormalized 
parameters will be contained in higher order corrections 
to the generalized interaction operator, and this correc- 
tions can be obtained order by order within the EFT ap- 
proach. This has been demonstrated in next-to-leading 
order by using the example of the EFT of short-range 
forces developed by van Kolck Q. 

We have shown that up to next-to-leading order the 
dynamics of the theory is governed by the generalized 
dynamical equation with the nonlocal-in-time interaction 
operator (74). These results have been obtained in the 
particular case where the parameters of the theory sat- 
isfy the condition (64). This limitation is a consequence 
of the fact that we used the assumption that the two- 
particle T-matrix of the theory is of the form (62), while 
this assumption does not fulfilled in general. The gener- 
alization of these results to the case where the condition 
(64) is not satisfied can only lead to another form of the 
generalized interaction operator. The general case will 
be discussed in detail in another paper. At the same 
time, the model we have considered in Sec.V provide 
the transparent illustration of the fact that in any or- 
der of the EFT approach low energy dynamics of nucle- 



ons is governed by the generalized dynamical equation 
(^) with a nonlocal-in-time interaction operator. As has 
been proved, the assumption that the relevant solution of 
Eq.(16) should be of the form (62) uniquely determines 
this solution and leads to the condition (64) . In this way 
we arrive at exactly the same on-shell T-matrix that has 
been obtained in Ref. by summing the bubble dia- 
grams up to next-to-leading order. The corresponding 
generalized interaction operator is of the form (74). The 
generalized dynamical equation (^ with this interaction 
operator is well defined and allows one to describe the 
dynamics of the system in a consistent way. 

The above method for constructing the effective inter- 
action operator is inapplicable in general, because the 
parameters of the theory need not satisfy the condition 
(64). In general, in order to construct the generalized 
interaction operator parametrizing the interaction in the 
EFT of nuclear forces, one has to analyze the contribu- 
tions from the off-shell amplitudes to T{z) in the limit 
|z| — > oo. In this limit the main contribution to the 
T-matrix comes from the bubble diagrams including po- 
tential pions that produce the pion-exchange potential. 
By summing the bubble diagrams in the large z limit, 
one can determine the contact component of the inter- 
action operator parametrizing the NN interaction. Thus 
the NN interaction operator can be represented as the 
sum (75) of the nonlocal contact and long-range com- 
ponents that can be obtained separately. The contact 
component is represented by a nonlocal-in-time interac- 
tion operator, while the long-range one is represented by 
the pion-exchange potential. In general, the long-range 
component should be supplemented by the Coulomb po- 
tential. It is extremely important, that starting with this 
NN interaction operator, we can construct the T-matrix 
and the evolution operator without resorting to summing 
all relevant diagrams and using regularization and renor- 
malization procedures. The generalized dynamical equa- 
tion (^ with this interaction operator is well defined and 
can be reduced to integral equations which do not re- 
quire renormalization. This has been demonstrated in 
leading order of the EFT approach. The operator (76) 
parametrizes the interaction that includes not only the 
leading order contact component but also the long-range 
one being described by the potential V{p2, Pi). The gen- 
eralized dynamical equation with the interaction opera- 
tor (76) can be reduced to the integral equation (78). 
This equation is so convenient for numerical calculations 
as the LS equation which is its particular case where the 
parameter Cq^"^ tends to zero and hence one can neglect 
the contact component. In this case the theory is reduced 
to the ordinary theory of the NN interaction based on the 
use of the pion-exchange potential. 

Finally, we have shown that the use of the formalism 
of the GQD allows one to formulate the effective field 
theory of nuclear forces as a completely consistent the- 
ory based on the well-defined equation of motion that 
does not require regularization and renormalization. Be- 
ing formulated in this way, the effective nuclear theory 
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permits a natural parametrization of the interaction of 
nucleons by the generahzed interaction operator that be 
derived order by order by using the methods of subtrac- 
tive eft's. One can hope that this operator will be able 
to play the same role in nuclear physics as the Coulomb 
potential in quantum mechanics of atomic phenomena. 
A remarkable feature of such a formulation of the theory 
is that in this case all advantages of the operator formal- 
ism of quantum mechanics can be used. The equation of 
motion (^), for example, allows one to construct not only 
the S-matrix but also the evolution operator describing 
the dynamics of nuclcon systems. This is very important, 
because the S-matrix is not everything. For example, at 
finite temperature there is no S-matrix because particles 
cannot get out to infinite distances from a collision with- 
out bumping into things [ p^ . In conclusion, it should 
be emphasized that the above is not a new approach to 
the EFT of nuclear forces. We have only shown that the 
EFT approach gives rise to low energy nuclcon dynamics 
which can be described, in a natural way, only by us- 
ing the generalized dynamical equation (||), and its use 
as the equation of motion describing low energy nuclcon 
dynamics can allow one to formulate the effective field 
theory of nuclear forces as a perfectly satisfactory theory 
like the quantum mechanics of atomic phenomena. 



APPENDIX A 

Let us consider the solution of Eq.(16) in the case 
where the interaction operator is of the form (76). From 
(17) and (18) it follows that this solution can be repre- 
sented in the form 



<P2|T(z)|pi >= lim <p2|r„(z)|pi >, (Al) 

u — > — oo 

where the operator Tu{z) is the sohition of the equation 
T„(z) = B{u) + (u - z)S(u)Go(u)Go(z)T„(2). (A2) 
Here the operator B{z) is given by 

<P2|B(z)|pi >-/l(z)+y(p2,Pl), 

with 

An 167r2 



(A3) 



/i(^) 



■'0 



The solution of Eq.(A2) can be represented in the form 

<P2|T„(z)|pi >-4")(z) + 4")(2;pi) 

+t<"'(z;P2)+4")(z;pi,P2). (A4) 

Substituting this representation in Eq. ( [A2| ) yields the 
following equations for t''Q^\z), t^j^\z;p),ti \z;p2) and 
4"^(2;pi,P2): 

C\z) = h{u) + {u-z)h{u) 

(4"^(.) + 4"^(^;k))^ 



(2^)3 {z-Ek){u-Ek) 



(A5) 



t^^\z;p) = {u- z) 



d^k (4"^(^) + 4"^(^;k)) 

(2^)3 iz-Ek){u-Ek) 

xT/(k,p); (A6) 



4"^(^;p) = {u- z)fi{u) 



pi"^(^;p)- 



+ 4^)(z;k,p) 



(27r)3 



{z-Ek){u-Eu) 



(A7) 



iu) / \ / \/ , F d k 

ty{z;p2,Pi) = V^{p2,Pi) + {u~z)J j^-p^ 



y(k,pi). (A8) 



(4")(z;P2)+4"'(z;p2,k) 
{z~Ek){u~Ek) 
It is not difficult to verify that 

4")(z;p)ee4'')(z;p). 

By solving the above set of equations, one can obtain 

the functions t'^\z), t^"\z,p) and 4"''('Z: P2, Pi) that in 
turn can be used for constructing the T-matrix. In fact, 
from (Al) and (A2) it follows that the T-matrix can be 
represented in the form ( |7^ ) where the functions to{z), 
ti{z,p) and i2(^,P2,Pi) are given by 



to{z) ^ lim 4"'(^)' tii^^P) = 1™ t'{"''{z,p), 

u — !■ — OO u — > — oo 

t2(2,P2,Pi)= lim 4"^(2:P2,Pl)- 



Ml 



u — * — oo 



Taking into account that 

(Pk u — z 

(2^)3 {z-Ek){u-Ek) 



my/m 
47r 



Eq.(A5) can be rewritten in the form 



47r 



16^2 \ 



47r 



m2 ^/—u 
raJm , 



+ 4"^w 



3 

7712 - 



/ 47r 
4^ /■ ,3,ti(2,k) 



m2 \/ —u 



(Fk- 



Ek 



o[\u\ 



-1/2X 



Letting m — > — oo in this equation and assuming that 
y(p2, pi) satisfies the ordinary requirements of quantum 
mechanics, one can easily get Eq.(|8l|). In the same way, 
from (IaS) and (^ one can derive Eqs.(f8|) and (|3). 
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